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The  Similarity  Between  Shapes  under  AfHne  Transformation 

Hong,  Jiawei  and  Tan,  Xiaonan 
Courant  Institute,  New  York  University 

Abstract 

The  recognition  of  similarity  between  shapes  under  affine  transformation  is 
complicated  because  we  have  to  consider  shapes  with  rotational,  translational,  scal- 
ing or  stretching  degrees  of  freedom.  We  simplify  this  problem  by  introducing  the 
canonical  form  of  a  shape.  Any  shape  can  be  transformed  by  an  affine  transforma- 
tion to  a  canonical  form.  The  canonical  forms  of  all  shapes  which  are  identical  under 
affine  transformation  can  coincide  through  rotations.  A  systematical  way  is  provided 
to  find  the  measure  function  which  can  indicate  the  difference  between  shapes  and 
keep  independent  of  affine  transformation.  Then  we  can  reduce  the  recognition  of 
the  similarity  between  shapes  under  affine  transformation  to  the  computation  of  the 
affine-independent  difference  between  their  canonical  forms,  i.e.,  finding  the  best 
matching  between  a  pair  of  canonical  forms  by  rotation,  which  is  a  one  dimensional 
search  problem. 

We  have  designed  a  linear  algorithm  for  transforming  a  shape  to  its  canonical 
form.  The  computation  of  the  difference  between  a  pair  of  canonical  forms  can  be 
handled  by  various  current  techniques  and  a  new  method  is  presented  in  this  paper. 
This  method  is  robust  in  the  sense  that  small  changes  in  the  shapes  will  only  cause 
small  changes  in  the  final  results.  Not  only  similar  shapes  can  be  identified  but  also 
the  quantity  of  difference  between  shapes  can  be  computed. 

1.  Introduction. 

The  similarity  between  shapes  is  playing  an  important  role  in  computer  vision. 
Considerable  attention  has  been  paid  to  the  problem  using  the  similarity  of  shape 
characteristics  in  object  recognition. [2][3][7]]8][9][10]  However,  there  are  far  fewer 
discussions  on  the  same  problem  under  an  affine  transformation. [6]  It  is  common 
knowledge  that  an  image  of  an  object  obtained  by  human  eyes  or  cameras  is  actually 
a  projective  transformation  of  the  object.  When  the  distance  between  the  viewer  and 


the  object  is  not  very  short  a  projective  transformation  will  be  very  approximate  to 
an  affine  transformation.  The  problem  how  to  recognize  whether  two  shapes  are 
identical  or  similar  under  an  affine  transformation  is  practically  meaningful  and  dif- 
ficuh. 

If  we  can  make  two  shapes  identical  by  applying  an  affine  transformation  to  one 
of  them,  then  we  say  that  these  two  shapes  are  identical  under  affme  transforma- 
tion. To  recognize  the  similarity  between  shapes  under  affine  transformation,  first 
of  all  ,  we  have  to  find  some  measure  which  can  indicate  the  "difference"  between 
shapes  and  is  independent  of  affine  transformation,  i.e.,  the  difference  indicated 
keeps  invariant  after  the  shapes  undergo  affine  transformations.  (Most  existing 
measure  functions  do  not  have  this  property).  If  this  kind  of  difference  between 
shapes  is  small,  then  we  say  the  two  shapes  are  similar  under  affine  transformation. 
The  smaller  the  difference,  the  more  similar  the  shapes. [1]  Moreover,  an  affine 
transformation  can  be  decomposed  into  a  translation  (two  parameters),  one  rotation 
(one  parameter)  and  two  stretchings  along  a  pair  of  orthogonal  axes  (three  parame- 
ters). Therefore  the  task  of  finding  the  minimum  difference  between  two  shapes 
under  affine  transformation  suffers  from  a  search  problem  in  a  six  dimensional 
space. 

In  this  paper  we  propose  an  efficient  method  for  recognizing  the  similarity 
between  shapes  under  affine  transformation.  (We  consider  two  dimensional 
objects).  The  key  point  is  to  define  and  compute  the  canonical  form  of  an  arbitrary 
shape.  Just  like  the  principle  in  algebra,  when  we  want  to  know  whether  two 
matrices  are  identical  or  not  we  can  transform  these  two  matrices  to  their  canonical 
forms  to  check  if  the  canonical  forms  are  identical.  Using  our  method  any  shape  can 
be  transformed  by  an  affine  transformation  to  a  canonical  form,  which  is  centered  at 
the  origin  point  of  the  coordinates  and  whose  moment  along  any  axis  passing 
through  the  center  equals  to  1.  If  two  shapes  are  identical  under  affine  transforma- 
tion then  their  canonical  forms  can  coincide  through  rotation. 

Based  on  the  canonical  form  of  a  shape  we  develop  a  systematical  way  to  find 
the  functions  which  can  measure  the  difference  between  shapes  and  which  are 
independent  of  affine  transformation. 

By  introducing  the  canonical  form,  the  recognition  of  similarity  of  shapes  under 
affine  transformation  can  be  done  by  comparing  their  canonical  forms.  Then  this 
problem  can  be  reduced  to  a  one  dimensional  search  problem:  to  find  the  best  match 
between  a  pair  of  canonical  forms  by  rotation. 


Because  we  do  not  use  higher  order  moments,  which  are  very  unsteady,  our 
method  is  robust  in  the  sense  that  small  changes  in  the  shapes  will  only  cause  small 
changes  in  the  final  results  obtained. 

This  paper  is  organized  into  6  sections.  In  the  next  section  we  define  dissimilar- 
ity functions  which  measure  the  "difference"  between  shapes.  In  section  3,  we  dis- 
cuss the  dissimilarity  functions  which  are  independent  of  some  transformations 
(translation,  rotation,  scaling  and  affine  transformations).  We  introduce  the  canoni- 
cal form  of  a  shape  under  affine  transformation  in  section  4.  In  section  5,  we  give  a 
linear  algorithm  to  compute  the  canonical  form,  we  discuss  how  to  evaluate  the  dis- 
similarity in  section  6. 


2.  Similarity  and  dissimilarity 

We  compare  shapes  by  their  features.  The  closeness  of  features  reflects  close- 
ness of  shapes.  To  recognize  the  similarity  between  shapes  we  have  to  define  some- 
thing which  can  measure  the  difference  between  shape  features. 

Definition  2.1    If  for  any  pair  of  shapes  iP,Q)  we  have  defined  a  non-negative 
number  d{P,Q)  satisfying 

(1)  d{P,P)  =  0 

(2)  d{P,Q)  =  d{Q,P) 

(3)  diP,R)^d(P,Q)+d(Q,R) 

then  we  say  d{P,Q)  is  the  dissimilarity  between  P  and  Q. 

According  to  the  definition,  two  identical  shapes  are  similar;  if  shape  P  is  simi- 
lar to  Q  then  shape  Q  is  similar  to  P;  if  P  is  similar  to  Q  and  Q  is  similar  to  R  then  P 
is  somehow  similar  to  R.  It  is  the  case  in  practice. 

Example     2.1      If     we     use     area{P)     to     represent     the     area     of    P.     then 
\area{P)—area(Q)  |  is  a  dissimilarity  function. 

Example  2.2  Hausdorff  distance.  First  define  the  distance  of  a  point  a  from  a  set  P 
by 

Hd(a,P)  =  Min{\a-b  \\b^P} 

Then  Hausdorff  distance  is 


Max{Hd(a,Q)\aiP}  +  Max{Hd{b,P)\biQ}. 

Hausdorff  distance  is  a  dissimilarity  function. 

Example  2.3  areaiP  +  Q),  where  +  is  the  exclusive  or  operator,  is  a  total  dissimi- 
larity function. 

Now  we  discuss  some  systematical  methods  for  finding  dissimilarity  functions. 

In  example  2.1,  the  dissimilarity  function  \area(P)—area(Q)  \  is  produced  by  a 
mapping  P -^area(P)^R,  where  R  is  the  set  of  all  real  numbers,  a  distance  space. 
More  generally,  assume  that  S  is  an  arbitrary  distance  space  and  /  is  an  arbitrary 
mapping  that  maps  a  shape  P  to  a  point  in  5,  then  we  can  define  a  function 

Di(P,Q)  =  d(f(P),fiQ)) 

through  the  distance  function  d  in  the  space  S. 

Theorem  2.1  D\{P,Q)  is  a  dissimilarity  function. 

Proof.  (1)  DiiP,P)=d(fiP),fiP))  =  0 

(2)  Di(P,Q)  =  d{f{P),fiQ))  =  d(f(Q),f(P))  =  Di(Q,P) 

(3)Di(P,R)=d(f(P),f(R)) 

^d(f(P)J(Q))  +  d(f(Q),f(R)) 

=Di(P,Q)+Di(Q,R) 

By  this  theorem,  we  can  produce  a  dissimilarity  function  from  any  feature  of  a 
shape  like  the  perimeter,  the  diameter,  the  logarithm  of  area and  so  on. 

Another  general  method  for  obtaining  dissimilarity  function  is  as  follows. 

Theorem  2.2  U  d{P,Q)  is  a  dissimilarity  function  and /is  a  concave  non-decreasing 
function  satisfying  that /(0)  =  0,  then  D2(P,Q)=f(d(P,Q))  is  a  dissimilarity  func- 
tion. 

Proof.  (1)  D2iP,P)=fidiP,P))=f  iO)  =  0 

(2)D2{P,Q)=fid(P.Q))=f{d{Q,P))=D2(Q.P) 
(3)D2iP,R)=f{d(P,R)) 

^f(d(P,Q)  +  d(Q,R))  (because  /is  non-decreasing) 


=fid{P,Q)+d(Q,P))+f(0) 
^fid{P,Q))+f(diQ,R))  (because /is  concave) 

=D2iP,Q)+D2{Q,R) 

We  can  have  a  sequence  of  "weaker"  and  "stronger"  dissimilarity  functions 
between  pairs  of  shapes,  the  stronger  requiring  that  more  features  of  the  shape  be 
similar  before  they  are  close.  By  introducing  a  dissimilarity  combination  function  we 
can  take  more  features  of  a  shape  into  account,  and  make  the  dissimilarity  function 
stronger. 

Definition  2.4  A  dissimilarity  combination  function  c{d\,d2)  is  a  function  satisfying 
that  if  di(P,Q),  d2(P,Q)  are  dissimilarity  functions  then  so  is 
c(diiP,Q),d2(P,Q)). 

Example  2.4  For  any  p  >l,c(rfi,rf2)  =  (^i^i''  +  ^2^5)^'''  is  a  dissimilarity  combina- 
tion function,  where  Xi  and  X2  are  used  to  weight  di{P,Q),  d2iP,Q),  respectively. 


3.  Invariant  Dissimilarity  and  Independent  Dissimilarity  Functions. 

The  similarity  between  shapes  should  be  at  least  independent  of  translation, 
scaling  and  rotation.  Therefore,  a  good  dissimilarity  function  should  reflect  this  pro- 
perty. 

Definition  3.1.  Let  G  be  a  transformation  group  on  a  plane.   A  dissimilarity  function 
d(P,Q)  is  G  invariant  if 

d(7P,7Q)  =  d(P,Q)  for  aWriG. 

A  dissimilarity  function  d{P,Q)  is  G  independent  if 

rf(T'P,T"e)=^(P,(2)  for  all  t',t"€G. 

The  examples  of  dissimilarity  functions  in  the  last  section  are  all  translation, 
rotation  invariant,  but  not  scaling  invariant.  Both  Hausdorff  distance  and  function 
areaiP  +Q)  are  neither  translation  independent  nor  rotation  independent.  All  of 
them  are  not  scaling  independent. 

As  mentioned  earlier,  we  are  interested  in  recognizing  the  shape  similarity 
under  affine  transformation.  Under  an  affine  transformation  the  shape  of  an  object 
can  undergo  an  arbitrary  scaling,  translation,  rotation  and  stretching  (a  compression 


is  regarded  as  a  stretching  with  factor  less  than  1)  along  a  pair  of  orthogonal  axes. 
The  situation  is  more  complicated.  All  examples  in  the  last  section  are  neither  affine 
independent  nor  affine  invariant.  It  seems  not  easy  to  find  affine  invariant  dissimi- 
larity functions.  It  is  even  harder  to  fmd  affmc  independent  dissimilarity  ones. 

Here  we  give  an  example  of  affine  invariant  dissimilarity  function 

Example  3.1  d(P,Q)  =  area(P  +  Q)/area(P\JQ). 

It  is  not  trivial  to  prove  that  this  function  qualifies  for  an  affine  invariant  dis- 
similarity function.  Notice  that  this  example  is  not  affine  independent. 

we  can  prove  that  if  J  is  a  G  invariant  dissimilarity  function,  then 

D3(P.Q)  =  Min{diP.7Q)} 

is  a  G  independent  dissimilarity  function.  However,  because  Min{diP,iQ)}  is  diffi- 

cult  to  compute,  this  method  is  impractical.  In  the  next  section  we  will  present  a 
better  approach. 


4.  Canonical  Form  under  Affine  Transformation 

Before  discussing  the  canonical  form  of  a  shape  under  affine  transformation,  let 
us  introduce  the  moment  curve  of  a  shape  first. 


Fig.  1 


As  shown  in  Fig.  1,  point  O  is  the  center  of  gravity  of  a  shape  P.  and  L  is  a  line 
passing  through   O.  Choose  coordinate  axes  such  that  the  origin  point  O  coincides 


with  the  center  of  gravity  of  the  shape  P  and  X-axis  coincides  with  line  L.  The 
moment  of  P  is  the  integration 

e=Sy'^dxdy. 
p 

We  take  a  point  M  on  OX  such  that  its  distance  from  the  origin  point  O  equals  to 
this  integration.  We  rotate  the  coordinate  axes  about  the  origin  point  O  and  denote 
the  angle  between  line  L  and  the  X-axis  by  6.  Then  0  changes  from  0  to  Zir.  By  trac- 
ing point  M  on  X-axis,  we  obtain  a  curve.  By  a  well  known  conclusion  from 
geometry,  the  curve  can  be  expressed  as 

€=a-t-fcCo5(2e-l-4)),  0<e<2'ir,  a>:fe>0. 

This  curve  is  called  the  moment  curve. 

Now  we  introduce  a  very  useful  concept:  the  canonical  form  of  shapes  under 
affine  transformation. 

Definition  4.1   A  shape  is  in  canonical  form,  if  it  satisfies  that 

(1)  The  center  of  gravity  is  at  the  origin  point. 

(2)  Its  moment  curve  is  a  unit  circle. 

Theorem  4.1  Any  shape  with  positive  area  can  be  transformed  to  a  canonical  form 
by  an  affine  transformation. 

Proof.  The  moment  curve  can  be  expressed  as  e  =a  +b  cos{2Q  +  ^),  0^6^27r, 
a>t^O.  When  6=— <t>/2  the  value  of  e  will  reach  its  maximum,  while  when 
8  =  17/2  — 4>/2  the  value  of  e  will  reach  its  minimum.  These  two  angles  are  orthogo- 
nal. By  two  stretches  along  two  orthogonal  axes  corresponding  to  these  two  angles, 
we  can  transform  P  into  another  shape  P'  such  that  the  moments  of  P'  along  these 
two  directions  (8  =  — 4>/2  and  e  =  ir/2-<t>/2)  are  equal  to  1.  Let  the  moment  curve 
of  P'  be  expressed  by  e=A  -1-B  Co.y(28  +  <J)).   Then  we  have 

\  =  A  +fiCo.y(-<t)  +  4>)=A  -1-5. 

\=A  -l-fiCo5(iT-<t>-l-<t))=A  -B. 

Therefore  A  =  l  and  B  =  0.  The  equation  becomes 

e  =  \. 

This  means  the  moment  curve  of  P'  is  a  unit  circle.  To  transform  a  shape  to  its 
canonical   form   we   need   to   translate   the   shape   such   that   its   center   of  gravity 


coincides  with  the  origin  point,  and  stretch  the  shape  along  a  special  pair  of  orthogo- 
nal axes,  therefore,  any  shape  can  be  transformed  to  its  canonical  form  by  an  affine 
transformation. 

Theorem  4.2  If  P  and  Q  are  identical  under  affine  transformation  and  C(P),C(Q) 
are  their  canonical  forms,  then  C  {P)  and  C{Q)  can  coincide  through  a  rotation 
about  the  origin  point,  i.e.,  for  any  affine  transformation  t  there  is  a  rotation  r 
about  the  origin  point  such  that  C  {7P)  =  rC(P). 

Proof.  Affine  transformation  form  a  group.  If  t  is  an  affine  transformation  then  so 
is  T~^;  if  Ti  and  T2  are  affine  transformations  then  so  is  T1T2.  C(P)  can  be 
transformed  to  P,  P  can  be  transformed  to  Q  and  Q  can  be  transformed  to  C{Q),  by 
affine  transformations.  Therefore  C (P)  can  be  transformed  to  C((2)  by  an  affine 
transformation  t,  i.e.,  C (P)  can  coincide  with  C (Q)  through  t.  Any  transformation 
can  be  decomposed  into  a  translation,  two  stretches  and  a  rotation,  therefore  we  can 
write  T  =  ToKiK2r  where  tq  is  a  translation,  ki,K2  are  two  stretches  along  two  mutu- 
ally orthogonal  axes  passing  through  the  origin  point  O  and  r  is  a  rotation  centered 
at  the  origin  point.  Since  C  (P)  and  C (Q)  have  the  same  center  of  gravity  (  the  ori- 
gin ),  the  translation  tq  must  be  identical,  to  =  ^,  where  e  is  the  identity  transforma- 
tion. Because,  if  to=!c,  the  to  would  translate  the  origin  point  O  to  a  non-origin. 
But  Ki,K2,r  can  only  transform  origin  to  origin,  and  non-origin  to  non-origin,  there- 
fore T  =  ToKiK2r  would  transform  the  origin  point  to  a  non-origin,  a  contradiction. 
Since  C(P)  and  CiQ)  have  the  same  moment  curve  (  the  unit  circle),  we  must  have 
Ki  =  K2  =  e.  Because  if  any  of  them  were  not  the  identical  transformation  e,  then  they 
would  change  the  moment  curve  from  a  unit  circle  to  some  thing  else,  which  is 
impossible.  Therefore  T  =  eeer  =  r  and  C  {P)  can  be  transformed  by  a  rotation  T  =  r 
toC(Q). 

In  the  last  section  we  mentioned  that  it  is  not  easy  to  find  an  affine  invariant  or 
independent  dissimilarity  function.  But  now  based  on  the  canonical  form,  we  can 
find  many  of  them. 

Theorem  4.3  If  d(P,Q)  is  a  rotation  invariant  dissimilarity,  then 
DdP,Q)=MIN{d(C{P),rC(Q))} 

r 

is  an  affine  invariant  dissimilarity.  The  minimum  is  among  all  rotations. 
Proof.  It  is  a  dissimilarity  function  because: 


i\)D4(PP)=Mm{diC(P),rC(P))}^d(C{P),eCiP))  =  0. 

r 

(2)  D4(P,Q)  =  Min{d(CiP),rC(Q))} 

r 

=Min{d{r~'^CiP),r~'^rC(Q))}    (<i  is  rotation  invariant) 

r 

=  Min{d{r-'^CiP),C(Q))}    ir-h  =  e) 

r 

=  Min[d(C (Q) ,r-'^C(P))}    (d  is  symmetric) 

r 

=Min{d(C(Q),rC(P))}     (the  set  {r~^}  is  the  same  as  {r}) 

r 

=D4(Q.P). 

(3)  D4(P./?)  =  A//>i{^(C(P).rC(;?))} 

r 

=  A/m{J(r"^C(P),r~VC((2))}    (^  is  rotation  invariant) 

r 
=  Mm{J(/--^C(P).C(C))}     (r-V  =  e) 

r 

=  A/m{^(C(P),r'r"C(/?))}     (the  set  {r'r"}  is  the  same  as  {r}) 
'■'.'■" 

<A//n{J(C(P),r'C((2))+J(r'C(e),r'r"C(^))}     (triangular  inequality) 
'■','•" 

=  Min{d(CiP),r'C{Q))+diC(Q),r"C(R))}     (J  is  rotation  invariant) 

r',r' ' 

=  Min{d(C(P),r'C(Q))}  +  Min{d(C(Q),r"C(R))} 

r  r" 

=  D^{P,Q)+DdQ,R). 
It  is  affine  invariant,  because: 

04(T'P,T"e)  =  A///i{J(C(T'P),rC(T"e))}    (definition) 

r 

=  Mz>2{rf(r'C(P).rr"C(G))}    (Theorem  4.2) 

r 

=A/m{J(r'C(P),rVC(Q))}     (the  set  {r]r"  is  the  same  as  r'{r}) 

T 

=Min{d(C(P),rC(Q))}    (^  is  rotation  invariant) 

r 

=D4(P.Q). 


All  examples  of  dissimilarity  functions  presented  in  this  paper  are  rotation 
invariant.  From  each  of  them  we  can  obtain  an  affine  independent  dissimilarity 
function. 

Another  general  way  to  get  affine  invariant  dissimilarity  functions  is  as  follows. 

Theorem  4.4  Assume  that  S  is  a  distance  space  with  distance  function  d,  and  that  / 
is  a  rotation  invariant  mapping  which  maps  the  canonical  form  C{P)  to/(C(P))€5. 
Then 

D5(P.e)  =  ^(/(C(P),/(C(G))) 

is  an  affine  invariant  dissimilarity  function. 

Proof.  It  is  easy  to  show  D5  is  a  dissimilarity  function.  It  is  affine  invariant, 
because: 

D5(T'P.T"C)  =  ^(/(C(T'P))./(C(T"e)))     (definition) 

=  J(/(r'C(/')),/(r"C(e)))     (Theorem  4.2) 

=^(/(C(P)),/(C(e)))     (/is  rotation  invariant) 

=  Ds{P,Q). 

Theorem  4.5  If  Ji  and  d2  are  affine  invariant  dissimilarity  functions,  then  so  is  any 
combination  C{d\,d2). 

Proof.  Obvious. 

By  introducing  the  canonical  form  we  developed  a  systematical  way  to  find 
affine  independent  dissimilarity  functions.  In  fact,  we  can  prove  that  any  affine 
independent  dissimilarity  function  can  be  defined  this  way.  This  is  the  first  reason 
why  we  discuss  the  affine  canonical  form.  The  second  reason  is  that  using  the 
canonical  form,  the  computation  of  the  dissimilarity  between  two  shapes  under 
affine  transformation  can  be  considerably  simplified,  as  shown  in  the  next  two  sec- 
tions. 


5.  Find  the  Canonical  Form 

In  this  section,  we  show  how  to  compute  the  canonical  form. 


The  first  step  in  finding  the  canonical  form  is  to  compute  the  center  of  gravity 
and  translate  the  shape  to  make  the  center  of  gravity  coincided  with  the  origin. 

Now  it  is  necessary  to  compute  the  parameters  (the  angle  of  axis  and  two  fac- 
tors) of  the  stretches  that  will  change  the  moment  curve  to  a  unit  circle. 

The  moment  curve  is  expressed  by  equation 

m=a  +  bCos(2e  +  ^). 

Let    us    evaluate    the    moments    at    6  =  0, 211/3,417/3,    (see    Fig.    1)    and    obtain 
m\,m2,m-i.  We  have 

m  1  =  a  +  bCos  (j) 

m2  =  a+fcCo5(2iT/3  +  <}>) 

m3  =  <3 +feCo5(4iT/3  +  <t)) 

Adding  them  together,  we  obtain 

mi  +  wi2  +  '"3  =  3a, 

a  =  —{mi+m2+m2,). 


Squaring  them, 


-^^2  =  ^2^-^.2^  =  .^     ■    1 


(m2-a)2  =  fc2co52(2^  +  4.)  =  -^+yCo.y(4j  +  24») 
(;n3-^)^  =  fc^Co52(4Y  +  <t))  =  ^+jCo:f(2y  +  24>) 


Adding  them  together  we  obtain 


(mi-a)2-l-(m2-a)^  +  (m3-a)2  =  |-fc2 


m\—a 

<|)= arc  Co5 — ; . 

b 


When  6  =  — <|)/2,  m  arrives  its  maximum  value  a  +  b.    When  6  =  ir/2-<t)/2,  m 
arrives   its   minimum   value   a-b.     Therefore   after   stretching   by    ^a+b    along 


direction  6  =  --^  and  by  ^a-h  along  direction  8  =  iT/2-<t>/2,  the  shape  should  be 

in  a  canonical  form. 

It  is  easy  to  see  that  the  transformation  of  a  shape  to  its  canonical  form  can  be 
finished  in  0(«)  time,  where  n  is  the  total  number  of  points  or  line  segments  in  the 
representation  of  the  shape. 


6.  Computing  the  Dissimilarity  between  Shapes 

Having  changed  two  shapes  into  their  canonical  forms,  the  recognition  of  the 
similarity  under  affine  transformation  becomes  easier.  If  two  shapes  are  identical 
under  affine  transformation,  their  canonical  forms  will  coincide  through  a  rotation. 
If  we  want  to  detect  the  closeness  of  the  shapes  under  affine  transformation,  we  can 
choose  some  affine  independent  dissimilarity  functions  to  measure  the  difference. 

We  have  already  pointed  out  that  the  most  general  affine  invariant  dissimilarity 
function    is   Di,{P,Q)  =  Mm{d{C{P),rC{Q)))    .which    is    induced    from   a   rotation 

T 

invariant  dissimilarity  function  d(P,Q).  To  compute  the  value  of  D^{P,Q),  we 
should  find  the  rotation  r  that  transforms  C(Q)  to  a  best  match  position  with  C(P). 
This  is  a  minimizing  problem  in  one  dimensional  space. 

Many  methods  including  various  heuristic  algorithms  can  be  used  for  the  prob- 
lem of  best  match  in  one  dimensional  space.  In  this  section,  we  propose  a  new  one. 
Let 

*(r)  =  J(C(P).rC(Q)) 
The  function  ^{r)  usually  satisfies  Lipsiz  condition: 

|<I'(/-l)-<I>(r2)l^Llri-r2|. 

with  a  coefficient  L  depending  on  the  dissimilarity  function. 

Given  an  error  tolerance  €,  a  straightforward  way  to  compute  d{C iP),rCiQ)) 
is  to  try  all  values  of  r  with  a  step  length  €  /L.  Then  we  need  —  steps. 

In  fact,  we  can  make  use  of  the  property  |^(ri)-<I>(r2)  |^Z,  |ri  — r2 1  still 
further.  As  illustrated  in  Fig.  2,  we  draw  line  S  of  slope  —L  through  point  (a,^(a)) 
and  a  line  T  through  point  (fe,<l>(fe)).  These  two  lines  intersect  at  point  (u,v).  The 
Lipsiz  condition  guarantees  that  if  a^x^b  then  <I>(x)^v.  It  is  easy  to  show  that 
v=(^(a)  +  ^(b))/2-Lib-a)/2,  therefore 


^(x)^i^{a)  +  ^(b))/2-Lib-a)/2 


'^(X) 


Based  on  this  fact,  we  will  try  to  skip  as  many  computations  as  we  can. 

Now  we  divide  our  computation  into  phases.  At  phase  i,  we  compute  <I>(r)  for 
all  values  of  2A:,ir.it,  =  1/2', 3/2', 5/2'' (2'-l)/2'. 

In  other  words. 

Phase  0:  Jto=l; 

Phase  1;  iti=l/2; 

Phase  2:  it2  =  l/4,3/4; 

Phase  3:  A3 = 1 /8,3/8,5/8.7/«. 

Assume  that  we  have  finished  phase  i  and  m  is  the  minimum  value  of  $  that  we 
have  ever  found.  By  the  reason  mentioned  above,  we  can  easily  see  that  if  an  inter- 
val [;72',(;  +  l)/2']  has  the  property  that 

(cI)(;72')  +  <l'(0"  +  l)/2''))/2-L/2'  +  ^>:/n 

Then  we  can  guarantee  that  no  points  in  this  internal  can  have  a  lower  value  than  m, 
therefore  we  can  remove  all  these  kinds  of  intervals. 

At  phase  I'  +  l,  of  course,  we  need  not  try  points  inside  these  intervals.  The 
procedure  goes  on  until  l/2'<€/L. 


This  method  can  be  improved  upon  by  incorporating  heuristics.  At  the  very 
beginning  we  use  heuristics  to  find  some  promising  initial  position  match,  i.e.,  make 
the  initial  value  of  <I>(r)  small.  This  will  lead  to  skipping  many  point  tests.  For 
example,  when  the  best  matching  happens,  the  farthest  point  of  C(P)  from  the 
center  of  gravity  will  match  the  forthest  point  of  C(Q).  Therefore  we  should  first 
rotate  C(P)  and  C  iQ)  so  that  the  farthest  points  coincide,  then  start  the  test. 
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